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Equilibrium distribution method vs. Drift method 
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The study of microsystems and the development of nanotechnologies require new techniques to 
measure piconewton and femtonewton forces at microscopic and nanoscopic scales. Amongst the 
challenges, there is the need to deal with the ineluctable thermal noise, which, in the typical exper- 
imental situation of a spatial diffusion gradient, causes a spurious drift. This leads to a correction 
term when forces are estimated from drift measurements [Phys. Rev. Lett. 104, 170602 (2010)]. 
Here, we provide a systematic study of such effect comparing the forces acting on various Brownian 
particles derived from equilibrium distribution and drift measurements. We discuss the physical 
origin of the correction term, its dependence on wall distance, particle radius, and its relation to 
the convention used to solve the respective stochastic integrals. Such correction term becomes more 
significant for smaller particles and is predicted to be in the order of several piconewtons for particles 
the size of a biomolecule. 

PACS numbers: 05.40.-a; 07.10.Pz; 



I. INTRODUCTION 



The precise measurement of small forces plays a cen- 
tral role in science and technology. Apart from the in- 
strumental challenge of measuring forces in the pico- and 
femtonewton range, the underlying concept of how forces 
are determined between macroscopic objects cannot sim- 
ply be scaled down to the micro- and nanoscale. This 
is easily understood by considering, e.g., a micron-sized 
object suspended in a liquid environment. Due to its 
Brownian motion resulting from collisions with the sol- 
vent's molecules, inertial effects become largely negligible 
and the trajectory will look rather different compared to 
macroscopic objects [1] [2]. Due to the irregularities of 
this type of motion, various possible mathematical de- 
scriptions for the trajectory of a microscopic object ex- 
ist, e.g., the stochastic differential equations suggested by 
Ito and Stratonovich. Under many conditions all these 
descriptions lead to the identical physical interpretation. 
This is no longer true when the diffusion coefficient of 
the suspended particle becomes position-dependent, as 
this typically occurs close to other particles or to a wall. 
In such cases, different mathematical descriptions of the 
Brownian motion are not identical and at most one cor- 
rectly describes the physical reality. Although the need 
for such corrections - usually referred to as spurious drift 
- has been realized several decades ago, e.g., for numerical 
simulations [3HS]; they are seldom applied when analyz- 
ing experimental data [3 [S] . While such corrections are 
often small, we have recently demonstrated that in the 
presence of a spatial dependence of the diffusion coef- 
ficient the negligence of the correct spurious drift term 
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FIG. 1: Main techniques of measuring forces at the micro- 
scopic and nanoscopic scales, classified according to their 
force resolution and the working conditions - surface/bulk 
- for which they are best suited: atomic force microscopy 
(AFM) [10], photonic force microscopy (PFM) [TTHT5] . and to- 
tal internal reflection microscopy (TIRM) [14H16] . 



can lead not only to the wrong amplitude but even the 
wrong sign of the forces acting on a micron-sized object 
suspended in a liquid environment, a situation often en- 
countered, e.g., in biophysical experiments [9]. 

In this article, we provide a systematic study of how 
the presence of a position-dependent diffusion coefficient 
alters the interpretation of the forces acting on a colloidal 
particle. We demonstrate that only when the correct spu- 
rious drift term is subtracted, the measured forces are in 
agreement with theoretical predictions. In particular, we 
show that the magnitude of the spurious forces increases 
when considering smaller particles as often employed in 
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biophysical experiments: the spurious forces acting on a 
biomolecule with a diamter of about 10 nm can be on the 
order of several piconewtons. This value is comparable to 
the elastic forces that are commonly found in experiments 
probing the biomechanics of single molecules [T7J [TH] . 

In Section II we discuss the physical origin of the math- 
ematical ambiguity in the description of Brownian mo- 
tion in a diffusion gradient and how this ambiguity may 
result in different formulae for the estimation of a force 
from experimental data. We also show that there are es- 
sentially two ways of measuring forces at the nanoscale. 
The first is based on the equilibrium distribution of a 
particle (or particles) subjected to an a priori unknown 
force. The second takes advantage of the fact that an 
applied force results in a drift of the particle. In Section 
III we provide some detailed information on the TIRM 
technique and the materials we employed. In Section IV, 
we report the results of the force measurements on par- 
ticles of various sizes and materials suspended in water 
close to a wall under the action of electrostatic forces and 
effective gravity, quantifying in each case the value of the 
spurious forces. In Appendix A, finally, we clarify how to 
describe correctly such phenomena using both stochastic 
differential equations and Fokker-Planck equations. 

II. THEORY: FORCE MEASUREMENTS IN A 
OVERDAMPED SYSTEM 

The forces acting on a microscopic object immersed in 
a fluid medium can be assessed either by studying the un- 
derlying potential or by studying their effect on the ob- 
ject's trajectory (for an overview of force measurement 
techniques and their force resolution, see Fig. [lj. The 
first approach - to which we shall refer as equilibrium 
distribution method - requires sampling of the equilib- 
rium distribution. Accordingly, it can be only applied 
under conditions where the investigated system is at or 
close to thermodynamic equilibrium with a heat bath. 
The second method - to which we shall refer as drift 
method - does not require the object to be at or even 
close to thermal equilibrium. Therefore, it can be applied 
also to systems which are intrinsically out-of-equilibrium, 
e.g., molecular machines, transport through pores, DNA 
stretching [T71 [T5] . The latter method, however, requires 
to detect the object trajectory with high sampling rates, 
which can be technologically more challenging, in partic- 
ular when combined with a high spatial resolution. 

A. Equilibrium distribution method 

A microscopic object in contact with a thermal heat 
bath at constant temperature does not come to rest, but 
keeps on jiggling around due to the presence of thermal 
agitation. When the particle is subjected to an external 
potential well U(z), this leads to a Boltzmann distribu- 
tion p(z) = exp( — where k B is the Boltzmann 
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FIG. 2: Schematic procedure of measuring a force F(z) by the 
equilibrium distribution method. From the measured equilib- 
rium probability distribution p(z) (a) of a particle one obtains 
the potential energy distribution U(z) — — lnp(z) (b), which 
then gives the force F(z) = ~-§^U{z) (Eq.[l]) (c). 

constant and T is the temperature of the heat bath. 

Accordingly, it is possible to sample the steady-state 
position probability distribution p{z) by measuring a 
large number of uncorrelated object positions (Fig.[2ja)). 
The equilibrium potential can be derived by U[z) = 
-k B T\n(p(z)) (Fig.^b)) and the force (Fig.gc)) by 

= dU(z) = k B T dp{z) 
{ ' dz p{z) dz ' 1 ' 

Due to the exponential dependence of the probabil- 
ity distribution on the potential depth, in typical exper- 
iments only potential minima within typically less than 
~5fcflT are explored by the particle. While Eq. ([I]) can 
be only applied to equilibrium conditions, it should be 
mentioned that that it has recently been demonstrated 
that for the specific case of noncquilibrium steady states 
(NESS) a similar relation, i.e. a generalized Boltzmann 
distribution, can be derived for the stationary particle 
equilibrium distribution and the conservative part of the 
potential [HHD]. 

B. Drift method 

Since for a microscopic body suspended in a liquid 
medium viscous forces prevail by several orders of mag- 
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FIG. 3: Schematic view of the propagator p(Az; Zo, At), 
which gives the distribution of the particle position incre- 
ments Az from its initial position zo after a time step At. 
The distribution is Gaussian for sufficiently small At with 
standard deviation y / 2D(zo)At (Eq. Q). The force can be 
estimated from the measured drift according to Eq. dEJ. 



nitude over incrtial effects, a constant force F applied to 
a microscopic particle results in a constant drift veloc- 
ity v = F/j, where 7 is the object's friction coefficient. 
Since v = Az/ 'At can be retrieved from the measured 
particle displacement Az within time At, the force can 
be measured accordingly as 



F 



7 



Az 
At' 



(2) 



For large forces this obviously leads to an univocal result. 
However, when the drift force amplitude is comparable 
to the effect of the thermal noise, the measured particle 
displacement Az and, thus, the drift force vary between 
identical experiments, leading to a statistical distribution 
of the measured values (Fig. [3]) 



F(z) = 7 



Az 3 -(*) 
At 



(3) 



where Azj(z) denotes the j-th experimental value of the 
particle's displacement after time At. 

Although Eq. ([3| is key to measure forces under 
nonequilibrium conditions, it is only valid in situations 
where the diffusion coefficient D — kT/j of the object to 
which the force is applied is constant. When D becomes 
position-dependent, i.e. 



D(z) = 



(Az {z) - {Az 3 {z))f 
2At 



(4) 
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FIG. 4: In the presence of a diffusion gradient (a), the propa- 
gator p(Az; zo, At) is different depending on where D is eval- 
uated, i.e. at Zi n i t i a i (b), at z m j d dic (c), or at z flna j (d). In the 
latter two cases this results in a shift of the position distribu- 
tion mean - the spurious drift. 



Eq. ([3| must be corrected by an additional term (9} [21 



F(z) = 7 (js) 



Azj (z) 
At 



a'y(z) 



dD{z) 
dz 



(5) 



to which we shall refer as spurious force and which may 
depend on the specific choice of a € [0, 1]. 

In the case of systems that are strongly coupled to a 
heat bath, thermodynamic consistency requires that a = 
1 |22) . In particular, this is true for a Brownian particle, 
which is the system we have experimentally investigated. 
More generally, other values of a might be possible when 
describing other stochastic processes [2"3"j . 



1. The physical origin of the spurious force 

A qualitative physical understanding of the correction 
term in Eq. ^ can be gained considering the effect of 
a diffusion gradient (Fig. |4|a)) on a Brownian particle 
initially localized at position initial (orange dashed line 
in Fig. [4]) at time to- In the simplest picture, the particle 
diffusion results in a dichotomic movement either to the 
left or to the right with the same probability; therefore 
after time At the particle is displaced by initial ±\/2~D At. 
In a more realistic picture, the particle final position has 
a continuous probability distribution. In both cases, as- 
suming D constant, the final particle position distribu- 
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tion p(Az; zq, At) is symmetric, such as the histogram of 
Fig.gb). 

In the presence of a diffusion gradient (Fig. [4|a)) the 
value of D is obviously different at the initial and final po- 
sition and, therefore, the evaluation of the displacement 
is not univocal. Assuming D = -D(zintiai) ~~ to which we 
shall refer as ltd convention - p(Az; zq, At) is symmet- 
ric, as in the constant diffusion case (Fig.[4jb)). However, 
it could be argued that D should be averaged over the 
particle displacement; assuming thus D = -D(z m iddie) ~ 
Stratonovich convention - p(Az; zq, At) becomes asym- 
metric (Fig. |4]Jc)), because the particle displaces further 
when moving towards increasing diffusion. Finally, as- 
suming D = -D(zfi na i) - isothermal, anti-Ito or backwards- 
Ito convention - p(Az; zq, At) becomes even more asym- 
metric (Fig. gd)). The spurious drift, and the related 
spurious force, account for such asymmetry. 

In more general terms, we might assume D = (1 — 
a)-D(zinitiai) + a-D(zfinai) and a e [0,1], with a = 
corresponding to the Ito convention, a = 0.5 to the 
Stratonovich convention, and a = 1 to the isothermal 
convention. The latter one, in particular, turns out to be 
the only thermodynamically consistent one for a system 
coupled to a heat bath [22] , 



III. METHODS AND MATERIALS 
A. Total internal reflection microscopy (TIRM) 

The trajectories of colloidal particles close to a wall 
are measured with TIRM P^HTB] and a scheme of a typ- 
ical setup is presented in Fig. [5ja). A p-polarized laser 
beam (A = 658 nm) is totally internally reflected at a 
glass-liquid interface generating an evanescent field de- 
caying into the liquid. A spherical colloidal particle in 
the vicinity of the interface scatters the evanescent light. 
The scattering intensity shows a marked dependence on 
the particle position I(z) = Iq exp(— j3z) where (3 is the 
inverse evanescent decay length and Iq — I(z = 0) 
(Fig. [5](b)). The scattering intensity timeseries I(t) is 
collected by an objective and recorded by a photomulti- 
plier (Fig.[5](c)). Inverting the position-intensity relation, 
we finally obtain the particles trajectory z(t) with a spa- 
tial resolution of a few nanometers (Fig. |5jd)). Since 
our system is in thermal equilibrium, the measured par- 
ticle trajectory allows us to calculate the forces via both 
the equilibrium distribution and the drift method, as de- 
scribed above. 

Due to the fractal nature of the Brownian motion, the 
value of / A At*| ~ ) which has to be calculated for the drift 
force (see Eq. fah) strongly depends on the chosen time 
interval At. On the one hand At should be as short 
as possible, on the other hand it cannot be made arbi- 
trarily small due to the finite experimental acquisition 
frequency. A reasonable tradeoff is to choose At so small 
that the spatial variation of the drift force during such 
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FIG. 5: Total internal reflection microscopy (TIRM). (a) 
Sketch of a typical TIRM geometry with a single colloidal 
particle in front of a planar transparent wall. The particle 
is illuminated by the evanescent field created by total inter- 
nal reflection of a laser beam at the glass-fluid interface while 
performing Brownian motion. The scattered intensity / is 
measured using a photomultiplier (PMT). The forces acting 
on the particle are due to gravity G e a and electrostatic in- 
teractions F e i and are indicated by green arrows. When the 
relationship between the scattered intensity and distance I(z) 
is known (b) the intensity time-series I(t) (c) can be converted 
into the trajectory z(t) (d). 



a time step is negligible. In order to meet these con- 
ditions in our experiment, we have stepwise reduced At 
until the probability distribution of Az became Gaus- 
sian, as expected when the above condition holds. Be- 
cause under our conditions (see next paragraph) the spa- 
tial gradient of the force acting on the particle is not 
constant, the above condition on At is expected to vary 
with the particle-wall distance. This is shown in Fig. [6] 
where we plotted the probability distribution of Az for 
At = 2 ms and At = 20 ms obtained at two particle- wall 
distances z — 220 nm and z = 400 nm which correspond 
to the distances with the largest and smallest spatial vari- 
ation of the force acting on the particle. As can be seen, 
for At = 2 ms, we always obtain a Gaussian distribu- 
tion while for At = 20 ms deviations from a Gaussian 
fit (dashed line) are observed. In all measurements pre- 
sented in the following we have chosen At = 2 ms. 



B. Gravitational and electrostatic force 

In order to compare the measured forces with theory, it 
is important to have full knowledge about the interaction 
mechanisms of a colloidal particle close to a wall. For an 
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FIG. 6: Influence of At. At the distance zo = 220 nm 
At = 2 ms leads to a Gaussian distribution of the position 
increments (a), but for At = 20ms a clear deviation from 
a Gaussian distribution appears. At z = 400 nm, where the 
spatial variation of the drift is much smaller, the distribu- 
tion is again Gaussian for At = 2 ms (c) and remains almost 
Gaussian also for At = 20 ms. 



electrically charged dielectric colloidal sphere suspended 
in a solvent, the interaction forces have been demon- 
strated to be described by [T4HI6] 
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FIG. 7: Comparison of measured (bullets) and calculated 
(line) normalized vertical diffusion coefficient D±/Dse for an 
R — 400 nm particle as a function of the particle- wall separa- 
tion z. 



vertical component of the diffusion coefficient is |26j 



D±(z) 



D 



SIi 



sinh (a(z)) Y^, 
1 



(7) 



= 1 (2n- 



7t(w+l) 



l)(2n+3) 

and 



where l(z) = 

2 sinh ((27t+l)q(z)) + (2w+l) sinh (2a(z)) 
4 sinh 2 ((n+0.5)a(z))-(2n+l) 2 sinh 2 (a(z)) 

a{z) — cosh" 1 (l + j|) . Dj_(z) is zero at the wall 
and monotonically increases with z approaching the 
bulk value at a distance of several particle radii away 
from the wall. When calculating D±(z) from the par- 
ticle trajectories of our TIRM measurements according 
to Eq. |4), indeed we find good agreement with the 
theoretica prediction (Fig. [7]). 



F(z) = Be KZ - G off . (6) 

The first term is due to double layer forces with k^ 1 
the Debye length and B a prefactor depending on the 
surface charge densities of the particle and the wall (see 
Tab. III). The second term describes the effective gravi- 
tational contributions G c s — |7ri? 3 (p p — p s )g with p p 
and p s the particle and solvent density and g the gravi- 
tational acceleration constant. Under our conditions the 
additional contribution of van der Waals forces can be ne- 
glected since they become relevant only at much shorter 
particle- wall distances [MJ [35] . 

C. Diffusion in bulk and diffusion gradient in front 
of a wall 

The Stokes-Einstein diffusion coefficient of a spheri- 
cal colloidal particle immersed in a solvent is -Dgg = 
ksT '/6nriR, where 77 is the shear viscosity of the liquid. 
Close to a wall the bulk diffusion coefficient decreases due 
to hydrodynamic interactions. From the solution of the 
creeping flow equations for a spherical particle moving 
near a wall assuming nonslip boundary conditions, the 



D. Sample preparation and parameters 

As sample cell we used a 2 mm thick cuvette com- 
prising two optical flats separated by a spacer of silicon 
rubber. The cell was filled with clean deionized water 
containing a very small number of colloidal particles and 
150 — 300 /imM NaCl salt to adjust the Debye screening 
length. In order to vary the spatial gradients of the dif- 
fusion coefficient we used particles with different radii R 
and densities since they sample different ranges of z/R: 
R = 400 nm (titanium oxide) , R — 655 nm (melamin) , 
and R = 1180nm (polystyrene). For further details see 
Tab.U 

IV. RESULTS AND DISCUSSION 

In Fig. [8] we show as closed symbols the force (and 
as inset the corresponding potential U(z)) acting on a 
R = 400 nm titanium oxide particle as obtained by the 
equilibrium distribution method (Eq. ([lj). Since the 
forces as determined from this method are unambigu- 
ous, they will be considered as the true forces acting 
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TABLE I: Sample parameters. 



R [nm] 


Material 


ps [g/cm 3 ] 


k 1 [nm] 


B [ P N] 


NaCl H 


400 


titanium-oxide 


2.54 


25 


68 


150 
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melamin 


1.51 


18 
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1.05 


18 
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FIG. 8: Force derived by the equilibrium distribution method 
(squares) for the 400 nm radius particle (see Tab. |l| and the- 
oretical expectation (line) according to Eq. The inset 
shows the corresponding potential U(z). 
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FIG. 9: Force derived by the drift method for the very same 
particle as in Fig. [8] The force measured according to Eq. ([5| 
with a — 1 (bullets) coincides with the theory (line), while 
with a = (circles) and a = 0.5 (triangles) there is a clear 
disagreement. 



on the particle. This is also supported by the fact that 
the experimental data are in quantitative agreement with 
Eq. (solid line), where G c g and k" 1 are taken from 
the experimentally known parameters (Tab. [I]) while the 
prefactor B has been treated as a fit parameter. The 
value B = 68 pN is in agreement with other TIRM mea- 
surements under similar conditions [Po] . 

The symbols in Fig. [9] correspond to the force-distance 
relation as obtained from the drift method using Eq. (|]) 
for a — (bullets), a — 0.5 (triangles) and a = 1 
(open squares) . Since the gradient of D vanishes far away 
from the surface, the force dependence on a is most pro- 
nounced close to the wall but weakens at larger z where 
the curves will eventually merge (only beyond the maxi- 
mum distance sampled by the particle) . The forces deter- 
mined with a = 1, i.e. the isothermal convention, show 
good agreement with the Eq. (|6| (solid line, the same as 
in Fig. [8]). We want to emphasize that all other choices 
of a, in particular negligence of the noise-induced correc- 
tion (a = 0), lead to significant differences. Not only the 
magnitude but also the sign of the forces obtained with 
a = 0, 0.5 disagree with the true forces as obtained in 



Fig-i 

Similar measurements were also performed with other 
particles which arc capable to sample even smaller 
particle-wall distance normalized by the particle radius, 
i.e. z/R (cf. Figs. 12 and 13), where the gradient of 
the diffusion coefficient becomes larger. Fig. [10] shows 



the results for the R — 655 nm particle and Fig. 11 for 
the R = 1080 nm particle. In both cases, the equilib- 
rium distribution measurement (orange squares) agrees 
with Eq. Q with B as the only fit parameter and with 
a = 1 (blue open squares) . As before, all other choices of 
a show no agreement, as exemplarily plotted for a = 
(blue bullets). 

In Fig. [l2j the experimentally determined spurious 
force obtained by the difference between the forces de- 
rived with the uncorrected drift method (Eq. ^ or 
Eq. (pi) with a = 0) and the equilibrium potential method 
(Eq. (|l])) is plotted for different particle sizes as a func- 
tion of z/R. The data are in good agreement with the 
solid lines representing the theoretical predictions for the 
spurious force a~f4-D(z) with a = 1. The spurious force 
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FIG. 10: Comparison of forces obtained for the 655 nm radius 
particle (see Tab.[l|. For the equilibrium distribution method 
(closed squares) and drift method with correction term (open 
squares) the results are in agreement with the theoretical ex- 
pectation (line) according to Eq. Q. Neglecting the correc- 
tion term leads to a clear deviation (bullets). 



FIG. 12: Distance dependence of the theoretically calculated 
spurious force 07 j^D with a = 1 for various particle radii 
R (solid lines). Experimental data is shown for R = 4 00 nm 
(circles, cf. Fig. [9|, R = 655 nm (squares, cf. Fig. 10 1, and 
R = 1180 nm (triangles, cf. Fig.[TT|. 
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FIG. 11: Comparison of forces obtained for the 1180 nm ra- 
dius particle (see Tab. fill. For the equilibrium distribution 
method (closed squares) and drift method with correction 
term (open squares) the results are in agreement with the 
theoretical expectation (line) according to Eq. Q. Neglect- 
ing the correction term leads to a clear deviation (bullets). 
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FIG. 13: Distance dependence of the spurious drift 4zD for 
various particle radii R (solid lines). Experimental data is 
shown for R = 400 nm (circles, cf. Fig. |9l|, R = 655 nm 
(squares, cf. Fig.|To|), and R = 1180 nm (triangles, cf. Fig.ITT}. 



depends only on the particle radius, but is independent 
of p s . It increases for particles nearer to the wall and 
for smaller particles due to the larger D, reaching val- 
ues in the order of several piconewtons for particles with 
R = 10 nm, the size of a macromolecule. 



Dividing the spurious forces (Fig. 12 1 by 7, it is pos- 
sible to derive the spurious drift a4-D\z), which is plot- 
ted in Fig. 13 for a = 1. The lines corresponding to 
the theoretical predictions and the symbols representing 
the experimental data for the particles with R — 400 nm 
(orange bullets), R — 655 nm (green squares) and R = 
1080 nm (blue triangles) show good agreement. Due to 
the increasing gradient of the diffusion coefficient when 
approaching a surface, the spurious drift increases for 
shorter particle-wall distances and for smaller particles. 
Differently form the spurious force, which diverges for 
small z, the spurious drift reaches a maximum at z = 0. 



V. CONCLUSIONS 

The experiments presented in this article clearly 
demonstrate that the presence of noise-induced drift has 
to be considered in force measurements based on the 
drift; neglecting it, this can lead to artifacts, which may 
even suggest the wrong sign of the force. Such spurious 
forces become more significant for smaller systems and 
may reach the piconewton range for objects the size of a 
biomolecule, i.e. about 10 nanometers. We stress, fur- 
thermore, that a constant diffusion can be assumed only 
for a particle far from any boundary. Such boundaries 
are naturally introduced by surfaces or by other parti- 
cles in suspension, a situation that is typically met in 
experiments. 



Appendix A: Mathematical description of the 
spurious drift: stochastic differential equations and 
Fokker-Planck equation 



mogorov extension theorem, for a thorough mathemati- 
cal treatment see, e.g., Ref. [29] ). For the Markov pro- 
cesses, such as processes describing motion of Brownian 
particles, finite-dimensional distributions are in turn de- 
termined by the initial distribution and by the transition 
probability density mentioned above. This shows, at an 
abstract level, equivalence of Kolmogorov's general ap- 
proach and Smoluchowski-Einstein transition probability 
description. 

While the Kolmogorov extension theorem restores the 
concept of an individual path to its position of central 
importance in the description of a diffusion process, it 
still leaves out another fundamental physical ingredient: 
infinitesimal evolution law. This is furnished by the the- 
ory of stochastic differential equations (SDE), which de- 
scribe time evolution of individual paths. This theory, 
developed in the 1940's by Ito [5T] (and to some extent 
also by Gikhman and by Stratonovitch) presents diffusion 
as a motion of a particle according to a random dynam- 
ical law. Solutions of SDE are Markov processes and, at 
least in principle, transition probabilities are determined 
by the equation and, as we have seen, they determine the 
statistics of the trajectories. 

The above remarks are purported to explain the sta- 
tus of the Fokker-Planck equation [5T] and its place in 
the diffusion theory. This partial differential equation 
(in mathematics called "Kolmogorov forward equation" ) 
has the transition density as its solution. It is thus an in- 
finitesimal approach to diffusion at the level of the transi- 
tion densities. As explained above, knowing the solution 
of the Fokker-Planck equation gives one full knowledge of 
the diffusion process. It is thus of utmost importance for 
a successful implementation of the Fokker-Planck equa- 
tion method, given an SDE, to be able to produce the 
correct form of the Fokker-Planck equation from its coef- 
ficients. Before we address this point, we have to discuss 
the so-called Ito-Stratonovitch dilemma, which is at the 
root of the existence of the spurious drift, and deals with 
the interpretation of SDE [315] . 



The motion of a Brownian particle can be described in 
several ways. In this Appendix we clarify the relations 
between the various approaches, putting them in a his- 
torical context and demonstrating how the spurious drift 
emerges. We refer to Ref. [2T] for further details. 

The first diffusion theory was developed by Smolu- 
chowski [57] and, independently, Einstein [25] at the 
beginning of the 20th century. The fundamental ob- 
ject of this theory was the transition probability den- 
sity, p Sit (initial, Zfinai) from position z ini tiai to position 
^finai between times s and t > s. The individual par- 
ticle trajectories played no major role in this descrip- 
tion of the diffusion process. In the 1930's Kolmogorov 
showed that defining a probability measure on the path 
space is equivalent to specifying all finite-dimensional 
distributions of a stochastic process, i.e. all joint dis- 
tributions of the random variables z(t\), 2(^2)1 ■ ■ • , z(tk), 
where t\ < ti < • • • < tk are some time instants (Kol- 



1. Stochastic differential equation approach 

We consider the one-dimensional, time-independent 
SDE 

dz = ^4 dt + j2D(z) dW(t), (Al) 
7(2) 

where W[t) denotes a Wiener process, i.e. a stochas- 
tic process, whose increments are stationary, indepen- 
dent and normally distributed with W(t) — W(s) hav- 
ing mean zero and variance equal \t — s\. Such equa- 
tion describes, for example, the behavior of a Brownian 
particle in the presence of a position-dependent diffusion 
coefficient D(z). While only one-dimensional SDE are 
discussed here, these considerations can be straightfor- 
wardly generalized to the multidimensional case. 
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TABLE II: Correct use of SDE and Fokker-Planck equations for systems coupled to a heat bath. 



a SDE Fokker-Planck equation 






dz = -EHdt+ ^D(z)dt+^2D{z) dW 




1/2 


dz = ~^dt+ \§- z D{z) dt + y/2D{z) dW 




1 


dz = -^dt+ s/2D(z) dW 





Eq. dA2j) should be interpreted as the integral equation 



z(T) 



dt - 



y/2D(z) dW{t). (A2) 



However, due to the irregularity of the Wiener process, 
as well as the solution z(t) of the SDE, the second in- 
tegral on the right-hand side has several possible inter- 
pretations. More precisely, it is defined as a limit of 

integral sums, Y,n=o a ( z ( t j))( w ( t ]+i) ~ W(*i))> where 
tj are points dividing the interval [0, T] into N equal 
subintervals and the intermediate points are defined by 
tj = (1 — a)tj + atj + \. A crucial point is that the value 
of the limit, i.e. of the stochastic integral, depends on 
the choice of a. Thus, in every applied problem, in ad- 
dition to the SDE, we have to know the value of a for 
the mathematical model of the studied phenomenon to be 
well-defined. Common choices are a = (the Ito conven- 
tion), a = 0.5 (the Stratonovitch convention) and a = 1 
(the isothermal, anti-Ito or backwards-Ito convention). 

The different choices of a are connected to each other 
by a precise mathematical relationship. Namely, the 
above SDE with a given choice of a is equivalent to the 
Ito equation (a = 0) 

dz= 1 ^dt + a dI ^dt+J2D{z)dW\ a=0 . (A3) 
j(z) dz 

That is, an equation with any choice of a can be rewritten 



equivalently as an Ito equation, at the cost of adding an 
additional term a. dD J^ dt. This term has been called 
"spurious drift" - a confusing name since it may suggest 
its nonphysical character, while, as we have seen, this 
additional drift term may be fully observable in a real 
physical situation [9]. 

Here, in a nutshell, we see the central problem ad- 
dressed in this paper: the choice of a leads in the Ito form 
of the equation to an extra drift term, which vanishes 
only when the diffusion coefficient is constant. When D 
is position-dependent, at most one of these values of a 
can correctly describe the system; such parameter may 
depend on the system under investigation. However, if 
the system is strongly coupled to a heat bath, as is the 
case for a Brownian particle, the spurious drift, and the 
associated spurious force, are maximal, i.e. a = 1 in 
Eq. (A3 1. Of course, as shown in Tab. |h| other conven- 
tions may also be used adjusting the weight of the spuri- 
ous drift term; remarkably, for the anti-Ito or isothermal 
convention the correction term vanishes |22J. 



2. Fokker-Planck approach 

Since the Fokker-Planck equation is deterministic, its 
solution involves no randomness and is thus uniquely de- 
termined. However, given a SDE, the prescription for 
writing the associated Fokker-Planck equation depends 
on the adopted convention, as is shown in Tab. |Hj 
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